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We have performed chain-growth simulations of minimalistic hybrid lattice models for polymers 
interacting with interfaces of attractive solid substrates in order to gain insights into the conformational 
transitions of the polymers in the adsorption process. Primarily focusing on the dependence of the 
conformational behavior on temperature and solubility we obtained pseudophase diagrams with a 
detailed structure of conformational subphases. In the study of hydrophobic-polar peptides in the 
vicinity of different types of substrates, we found a noticeable substrate specificity of the assembly of 
hydrophobic domains in the conformations dominating the adsorption subphases. 
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INTRODUCTION 

Understanding molecular self-assembly at organic-inorganic interfaces is essential for the 
design of related future nanotechnological applications, e.g., microscopic sensory devices in 
biomedicine and nanoelectronic circuits. Recently, the enormous progress in the development 
of high-resolution equipment allowed experiments which revealed quite interesting properties 
of such hybrid interfaces as, e.g., the specific dependence of peptide adhesion to the type of 
attractive substrates and peptide sequences [HEllS]. In studies of short peptides consisting 
of 12 amino acids, it was found, e.g., that the adhesion strength to a (100) siHcon (Si) 
surface improved by a factor of about 15 only by permuting the order of amino acids in 
this sequence. On the other hand, the adsorption strengths of the same sequences to gallium 
arsenide (GaAs) with (100) orientation hardly differ [2]. The reasons for this binding specificity 
are not yet understood and an appropriate atomic model explaining the specific substrate- 
peptide cooperativity on the microscopic scale is still lacking. This problem is related to similar 
studies where the adsorption and docking behavior of polymers is essential, e.g., protein-ligand 
binding [4], prewetting and layering transitions in polymer solutions as well as dewetting of 
polymer films |5j, molecular pattern, electrophoretic polymer deposition and growth [6]. 
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Fig. 1: Contact density g{ns,nm) of a homopolymer with 179 monomers in a cavity with 
attractive substrate. For regularizing the number of unbound conformations (i.e., rtg = 0), a 
steric wall parallel to the attractive substrate was placed in a distance of = 200 lattice 
units. 



CONFORMATIONAL TRANSITIONS ACCOMPANYING HOMOPOLYMER 

ADHESION 

For the following study of the conformational behavior of a homopolymer in the adsorption 
process to a solid, attractive substrate, we employ a strongly simplified coarse-grained model. 
The homopolymer is modeled as interacting self-avoiding walk on a simple cubic (sc) lattice 
in impHcit solvent, and the energy of a given conformation is related to the number Um of 
nearest-neighbor contacts of monomers being nonadjacent along the chain. In the adsorption 
phase, the surface contact energy in our model is proportional to the number of monomers in 
the surface-contacting layer, and the number of monomer-substrate contacts is denoted as n^. 
Introducing an overall energy scale eo (which is set to unity in the following) , and a parameter 
s that rates the energy scales of the polymer conformation and the surface contact energy, 
the minimaHstic model is written as [7l[8l[9]: 

Es{ns,nm) = -eQ{ns + sum)- (1) 

Since the parameter s effectively controls the compactness of the polymer conformations, it 
can also be interpreted as kind of solubility (the larger the value of s, the worse the quality of 
the solvent). In our simulations, we applied the contact-density chain-growth method which is 
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Fig. 2: Profile of the specific heat as a function of temperature T and solubility s. Black and 
white lines emphasize pseudophase transitions (cf. text). The dashed black Hne indicates the 
positions, where the free-energy minima of adsorbed and desorbed phases are degenerate, i.e., 
where both phases coexist with equal probability density. 

a generalized variant of the multicanonical chain- growth algorithm [lO]. This method allows 
a precise estimation of the contact density g{ns,nm), which is the number of conformations 
with Us surface and rim intrinsic contacts. For the homopolymer with 179 monomers used in 
our study, g{ns,n„i) ranges over more than 120 orders of magnitude, as shown in Fig.[TJ 

The main advantage of directly sampling g{ns,rLm) in our simulation is that this quantity 
is independent of the external parameters temperature T and solubility s. This means, the 
calculation of energetic statistical quantities can be performed after the simulation for arbitrary 
values of T and s. Defining the partition sum hy Z = ^ g(ns, rim) exp(— ii^s/fesT), 
ensemble averages of functions of the contact numbers f{ns,nm) are obtained via {f{ns,nm)) = 

n /("-si nm)g{ns, n-m) exp(— iJ^/fesT) /Z. As our main interest is focused on the conformational 
transitions experienced by the polymer in the adsorption process, we consider in the following 
the specific heat, which is here a function of temperature and solubility: Cv(T, s) — {{E"^) — 
{Es)'^)/kBT'^ ■ The profile of Cy is shown in Fig. [21 where bright regions high-light strong 
fluctuations, and black and white lines indicate conformational transitions. Black lines mark 
transitions that are assumed to resist the thermodynamic limit, whereas subphases speciflc to 
the precise length of the polymer are separated by white lines. The Cy profile can therefore 
also be interpreted as T-s pseudophase diagram of the hybrid system. There are two main 
regions, the desorption phases of respective compact and expanded conformations, DC and 
DE, and the phases of adsorbed conformations. The adsorption regime can be divided into 
compact film-like (ACl) and layered (AC2a-d) phases, globular conformations with surface 
contact (AG, AGe), and extended conformations (AE). In the low-temperature and poor- 
solvent pseudophases ACl, AC2a-d and AGe, layered conformations dominate. In ACl a two- 
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dimensional, maximally compact polymer film is entirely in contact with the substrate. The 
less the influence of the solvent, the higher is the tendency of the polymer to form a maximum 
number of intrinsic contacts by forming layers. In AC2d and AGe, double-layer conformations 
dominate, in AC2c triple-layer, and in AC2b four-layer structures. The maximum number 
of intrinsic contacts in three dimensions is found in the five-layer structures in subphases 
AC2ai,2- Note that trivial cubic symmetries are impossible as 179 is a prime number [UlS]. 

SUBSTRATE SPECIFICITY OF PEPTIDE ADSORPTION 

Considering peptides in the simplest representation, i.e., as a heteropolymer sequence of 
hydrophobic and polar residues on a simple cubic lattice, model ([1]), with slight modifications, 
can also be used to study peptide adhesion to substrates. In tertiary protein folding, the 
hydrophobic monomers form a compact core due to the hydrophobic effect, and the polar 
monomers screen this core from the aqueous environment. The two types of monomers naturally 
allow the investigation of three types of substrates: (a) the unspecifically attractive, (b) the 
hydrophobic, and (c) the polar substrate, where the latter two are only attractive to respective 
hydrophobic and polar monomers with surface contact. Therefore, model |T|) is substituted 
by [n]: 

i?s("-s, "-hh) = -£o("-s + stihh), (2) 

where tt-hh is the number of intrinsic nearest-neighbor contacts between nonadjacent hydrophobic 
monomers and Ug symbolizes the number of attractive surface contacts and thus depends 
on the type of substrate. In contact-density chain-growth simulations, we have studied the 
adsorption behavior of an exemplified heteropolymer of 103 monomers (66 polar, 37 hydrophobic), 
which is the hydrophobic-polar transcription of the amino acid sequence of cytochrome c. In 
Fig. [21 the specific heat profiles of the heteropolymer in the vicinity of the three different 
substrates are shown. Ridges (marked by white and gray lines) indicate conformational pseudophase 
transitions. In the bulk phases, the typical expanded random-coil-like conformations (DE) 
and the compact, native-like folds (DC) can be distinguished. In the adsorbed regime we also 
find expanded (AE) and compact/globular (AC, AG) phases. Typical conformations in the 
fine-structured AC subphases are also shown in Fig. [31 The compactness of the hydrophobic 
domains in these subphases does not only depend on the solvent quality, but also on the 
effect how polar residues hinder the formation of hydrophobic domains. On the hydrophobic 
substrate, there is an effective steric repulsion of the polar monomers which are pushed off 
the surface layer. In the case of active attraction of polar residues to the polar substrate, the 
competing tendency of the hydrophobic monomers to form compact, layered clusters in poor 
solvent leads to the loss of surface contact of these cores. 

An alternative glance at the phase behavior is thrown by exploring the free-energy landscape 
of the system. Considering the numbers of hydrophobic contacts nnn and surface contacts 
Us as natural system parameters, the free energy, expressed as function of these quantities, is 
defined as fr,s("-s, ^^hh) = -E's(ns,nHH) — TS{ns,nYi}i), where the "microcontact" entropy is 
related to the contact density (?(ns,7iHH) via S'(ns,nHH) — fcs In 5(71.^, tihh)- The minimum 
of Ft,s (f^s, ?^hh) for given external parameters s and T is related to a class of macrostates 
with nf^ surface and hydrophobic contacts which dominates the phase. In Fig. [H we 
have plotted the map of all possible free-energy minima in the range of external parameters 
T g [0, 10] and s G [—2, 10] for the 103-mer in the vicinity of an unspecifically attractive 
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Fig. 3: Specific-heat profiles for substrates being (a) unspecifically attractive, (b) hydrophobic, 
and (c) polar. Also shown are typical conformations in the AC subphases. 
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Fig. 4: Contact-number map of all free-energy minima for the 103-mer and substrate equally 
attractive to all monomers in the parameter space T e [0,10], s £ [—2,10]. Lines illustrate 
contact free energy changes with the temperature at constant solvent parameter s. For the 
exemplified solvent with s = 2.5, probability distributions close to the layering transition 
(AC1,2) and the unbinding transition (AE to DE) are also shown. 



substrate. Solid lines connect minima in the free-energy landscape when changing temperature 
under constant solvent (s = const) conditions. Following the exemplified trajectory for s = 2.5 
and starting at very low temperatures, it is clear from Fig. [SKa) that the system resides in 
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pseudophase ACl, i.e., compact, film-like single-layer conformations dominate. The system 
obviously prefers surface contacts at the expense of hydrophobic contacts. Increasing the 
temperature, the system experiences close to T « 0.35 a first-order-like conformational 
transition, and a second layer forms (AC2). The loss of energetically favored substrate contacts 
of polar monomers is partly compensated by the energetic gain due to the more compact 
hydrophobic domains. Increasing the temperature further, globular, pancake- like conformations 
dominate in the globular pseudophase AG. Reaching AE, the number of hydrophobic contacts 
decreases further. Extended, dissolved conformations dominate. The transitions from AC2 
to AE via AG are comparatively "smooth" (second-order-like), i.e., no immediate changes 
in the contact numbers passing the transition lines are noticed. The situation is different 
when approaching the unbinding transition line from AE close to T w 2.14. This transition is 
accompanied by a dramatic loss of substrate contacts - the peptide desorbs from the substrate. 
As the probability distribution in Fig. |4] shows, the unbinding transition looks again first- 
order- like, i.e., close to the transition line, there is a coexistence of adsorbed and desorbed 
conformations. 

CONCLUDING REMARKS 

We have studied the conformational behavior of minimalistic hybrid interfaces of polymers and 
substrates and obtained from sophisticated chain-growth simulations pseudo-phase diagrams 
which exhibit a rich pseudophase structure, in particular in the adsorption regime. We could 
also show that the adsorption of heteropolymers is specifically dependent on the hydrophobic 
or polar character of the substrate. These results, whose experimental verification is still 
pending, are of particular interest for future applications of hybrid organic-inorganic materials 
in nanotechnology and biomedicine. 
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